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ABSTRACT 

We use the spectra of «22,000 nearby early-type galaxies from the Sloan Digital Sky Survey (SDSS) to 
determine the age distribution of these galaxies as a function of their velocity dispersion a v in the range 
100 km s" 1 <, a v £ 280 km s . We then combine the inferred age-distributions with the local abundance of 
spheroids, including early-type galaxies and late-type bulges, to predict the evolution of the quasar luminosity 
function (LF) in the redshift range < z ^ 6. We make the following simple assumptions: (i) the formation of stars 
in each galaxy, at the epoch identified with the mean mass-weighted stellar age, is accompanied by the prompt as- 
sembly of the nuclear supermassive black hole (SMBH); (ii) the mass of the SMBH obeys the M^-Uy correlation 
observed in nearby galaxies; (iii) the SMBH radiates at a fraction /edd of the Eddington limit for a fixed duration 
/q, and is identified as a luminous quasar during this epoch, (iv) the intrinsic dispersions in the Eddington ratio 
and the M\,h-a v relation produce a combined scatter of AlogL<j around the mean logarithmic quasar luminosity 
(logLq) at fixed a v . These assumptions require that the SMBH remnants of quasars with L\, \ <, 10 12,5 /EddL© 
reside predominantly in bulges of late type galaxies. We find that evolution of the observed quasar LF can 
be fit over the entire redshift range in this simple model, < z ^ 6 with the choices of 0.6 ^ AlogLQ <S 0.9, 
6 x 10 7 yr <, ?q <, 8 x 10 7 yr, and 0.3 (/Edd) ^ 0.5. We find no evidence that any of the model parameters evolves 
with redshift, supporting the strong connection between the formation of stars and nuclear SMBHs in spheroids. 

Subject headings: quasars: general - galaxies: nuclei - galaxies: active - black hole physics - accretion 



1. INTRODUCTION 

The discovery of tight correlations between the masses of su- 
permassive black holes (SMBHs) at the centers of galaxies and 
the global properties of the spheroid component of the galax- 
ies themselves (see, e.g., Magorrian et al. 1998; Ferrarese & 
Merritt 2000; Gebhardt et al. 2000; Graham et al. 2001) have 
suggested a strong link between the formation of SMBHs and 
their host spheroids. Several groups have noted that galaxy for- 
mation and SMBH growth should be linked, and many mod- 
eled the joint cosmological evolution of quasars and galaxies 
(see, e.g., Monaco et al. 2000; Kauffmann & Haehnelt 2001; 
Haiman & Menou 2001; Granato et al. 2001; Ciotti & van Al- 
bada 2001, Cavaliere & Vittorini 2002; Lapi et al. 2006, and 
references therein). 

A particularly revealing characteristic of this cosmic evolu- 
tion is the so-called "downsizing". Observations have long 
suggested an inverse correlation between the ages of ellipti- 
cal galaxies and their size, with more massive ellipticals form- 
ing at earlier epochs (e.g., Heavens et al. 2004 and references 
therein). Interestingly, the characteristic quasar luminosity sim- 
ilarly declines over time from z ~ 2 to z = (Boyle et al. 2000). 
At first sight, the decline of both the characteristic galaxy veloc- 
ity dispersion, and of quasar luminosity, is surprising in "bottom- 
up" hierarchical structure formation models. This anti-hierarchical 
behavior, can nevertheless be naturally accommodated in these 
structure formation models, invoking the energy feedback from 
the AGN on the host galaxy (e.g. Granato et al. 2004) and 
large-scale environmental effects (e.g. Blanton et al. 1999). 
The similarity between the "downsizing" observed for elliptical 
galaxies and for quasars has been noted in previous work (e.g. 
Franceschini et al. 1999; Haiman, Ciotti & Ostriker 2004), and, 
as these works pointed out, it suggests that the formation of 



SMBHs and spheroids was linked by the same physical mech- 
anism at all cosmic epochs. 

The simplest hypothesis to interpret the above trends is that 
star formation in spheroids and BH fueling are proportional to 
one another at all times (Haiman, Ciotti & Ostriker 2004). With 
the additional assumptions that BH accretion luminosity stays 
at some fraction of the Eddington limit during luminous quasar 
phases, and that these phases have a characteristic duty cycle, 
one can predict the evolution of quasars from the evolution of 
early-type galaxies (and vice versa). Unfortunately, a direct 
observational test of the above hypotheses is not possible, as 
the cosmic evolution of the population of early-type galaxies 
has not been determined from observations out to high z- 

Here we follow an alternative "archaeological" approach, and 
we use age-distributions we infer for nearby early-type galax- 
ies, as a function of their velocity dispersion, to reconstruct in 
detail the evolution of quasars, as a function of their luminos- 
ity. This approach has been attempted previously by Cattaneo 
& Bernardi (2003), who have used the mean age of early-type 
galaxies us their velocity dispersion, together with assumptions 
about the geometry of obscuring tori, to infer the evolution 
of their abundance, and to predict the evolution of the quasar 
LF (see also Granato et al. 2001 for the inverse approach). 
In this paper, we improve on these previous studies in sev- 
eral ways: (i) we infer the age distribution, rather than only the 
mean age, as a function of er,,, for early type galaxies, (ii) we in- 
clude the population of late-type bulges in our analysis, which 
indeed dominate the spheroid population at low velocity dis- 
persion, (iii) we perform a quantitative statistical comparison 
between the predicted and observed quasar LFs, and (iii) the 
quasar lifetime, the mean Eddington ratio of SMBHs, and its 
scatter, are taken as three fitting parameters that are fit simulta- 
neously. We find that the shape of the quasar LF we predict is in 
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rough agreement with observations, without introducing either 
a luminosity-dependent obscuration or a redshift-evolution of 
either of the three model parameters. Our results require that 
the quasar luminosity Lq has a non-negligible scatter (of J> 0.5 
dex) at a fixed value a v of the host galaxy. 

The rest of this paper is organized as follows. In § 2, we dis- 
cuss the sample of early type galaxies we used in our analysis, 
and infer the age distribution of these galaxies. In § 3, we dis- 
cuss our modeling that converts the inferred age-distributions 
to the evolution of the quasar luminosity function (LF). In § 4, 
we discuss our method of statistically comparing the predicted 
and observed quasar LFs. In § 5, we present and discuss our 
main results, and in § 6, we briefly summarize our conclusions. 
Throughout this paper, we adopt the background cosmologi- 
cal parameters Q m = 0.3, Qa = 0.7, and Hq = 70 km s" 1 Mpc -1 , 
consistent with the values measured recently by the WMAP ex- 
periment (Spergel et al. 2006). 

2. THE AGE DISTRIBUTION OF EARLY-TYPE GALAXIES 

We have used the sample of early-type galaxies obtained 
by Bernardi et al. (2006a). The sample, extracted from the 
Sloan Digital Sky Survey (York et al. 2000), contains over 
40,000 early-type galaxies selected for having apparent magni- 
tude 14.5 < g < 17.75 with spectroscopic parameter eclass < 0, 
which gives a PCA component corresponding to no emission 
lines, and fracDev r > 0.8, which is a seeing-corrected indica- 
tor of morphology. Only those objects were included that have 
measured velocity dispersions, which translates in the spectra 
having S/N > 10 in the region 4200-5800 A. The sample ex- 
tends over a redshift range 0.013 < z < 0.25, which corresponds 
to a maximum look-back time of 3 Gyr. 

We computed the ages of the galaxies in three different ways. 
First, for each individual galaxy, we fitted a single stellar popu- 
lation model to the whole observed spectrum (continuum+lines) 
with two parameters: age and metallicity. Second, we used 
the MOPED algorithm (Heavens et al. 2000) to determine the 
star formation history of the galaxies. Third, we used the pub- 
lished ages by Bernardi et al. (2006a) that use stacked spectra of 
galaxies with similar properties to measure the absorption fea- 
ture produced by the Balmer absorption lines (e.g. H-f F ). Ages 
are then computed from the lines assuming a-enhancement 
with respect to the solar abundance stellar population models. 
The first and third method return a luminosity-weighted age, 
while the MOPED technique, because it reconstructs the star 
formation history, returns the age of the different stellar com- 
ponents, which can be used to computed the mass- weighted 
age for the galaxy as a whole. The three methods are compared 
in detail in a companion paper by Jimenez et al. (2006). 

There are about 22,000 galaxies in common between the 
MOPED sample (Panter et al. 2006) and the Bernardi et al. 
(2006a) sample. In Fig. 1, we show the age distribution for 
several velocity dispersions that we obtained for these galax- 
ies. The solid (black) histograms show the mass-weighted age 
distributions from MOPED. The ages correspond to those ex- 
trapolated to z = 0; note that we impose a prior that the age of 
any galaxy, when extrapolated to z = 0, cannot exceed 13.7 Gyr 
(Spergel et al. 2006). The distributions are broad for all values 
of a v , but are narrower for higher a v and centered on older ages. 
The uncertainties due to age-measurement errors are shown by 
the error bars in the figure, and will be discussed in § 4 below. 

For comparison, the blue (dashed) histograms show the age- 
distributions obtained from absorption line indices. As the fig- 
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FIG. 1. — The distribution of stellar ages of early-type galaxies, as a func- 
tion of their velocity dispersion a v in the range 100 km s" 1 < a v < 280 km s _1 . 
The spectra of ~ 22,000 galaxies, all from z < 0.25, were used from SDSS 
DR4. Each <T„-bin contains between ~300 and ~ 4500 galaxies, as labeled. 
The solid (black) histograms show the mass-weighted distributions inferred 
from the full spectra (continuum+lines) that we utilized in our analysis. Note 
that galaxies with larger <x v are not only older, they have a narrower age- 
distribution. The dashed (blue) histograms show, for comparison, the age- 
distributions derived from absorption line indices. This method generally gives 
shorter ages (see text for discussion). 



ure shows, the mean ages are lower, and the shapes of the dis- 
tributions are different, from the MOPED results. These dif- 
ferences are partly accounted for by the mass vs. luminosity 
weighting, and by the fact that the MOPED ages do not in- 
clude a correction for a-enhancement. Full spectral models 
for the continuum+line emission, including a-enhancement, 
are not available, but we expect that such models would some- 
what reduce the inferred ages. There appear to be additional 
differences between the two age-measurement techniques (see 
Jimenez et al. 2006 for further discussion). We chose to adopt 
the MOPED age-distributions in our analysis below, because 
they utilize more information from the spectrum, and because 
they can be used to compute mass-weighted ages. We will dis- 
cuss in § 5 below how our results would change if the ages 
inferred from the absorption-line indices were used instead. 

3. RECONSTRUCTING THE QUASAR LUMINOSITY FUNCTION 

We use the age-distribution of the early-type galaxies, in- 
ferred in § 2 above, to reconstruct the evolution of the quasar 
luminosity function, using the following simple assumptions. 

We start with the velocity function of early-type galaxies at 
z « 0, as determined by Sheth et al. (2003), 
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This function, shown by the dashed curve in Figure 2, shows 
a clear peak: the abundance of early type galaxies declines at 
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FIG. 2. — The local velocity function spheroids, adapted from Sheth et al. 
(2003), showing the contribution of early type galaxies (dashed curve), and 
the whole spheroid population, including the bulges of late type galaxies (solid 
curve). For comparison, the data-points show the quasar luminosity function 
at redshift z = 2.18 (adapted from Hopkins et al. 2006). The upper horizontal 
axis labels show the mean quasar luminosity corresponding to a v from equa- 
tion (3), which assumes that the BHs follow the M bh -o\, relation, and shine at 
a fraction /gad of the Eddington luminosity. 



velocity dispersions below <, 200 km s" 1 . Below this value, 
the population of spheroids is dominated by the bulges of late- 
type galaxies, as demonstrated by the solid curve showing the 
sum of the two spheroid velocity functions. For reference, the 
data-points in the figure show the observed quasar LF at red- 
shift z = 2.18 (adapted from Hopkins et al. 2006), using the 
mean Lq - a v relation from equation (3) without any scatter. It 
is clear from the figure that the quasar LF at faint luminosities is 
monotonic, and traces the total spheroid population much more 
closely than the subset represented by the early type galaxies. 
The empirical evidence listed in the Introduction also suggests 
that the correlation between SMBH mass and spheroid proper- 
ties apply both to SMBHs in early type galaxies and to those in 
bulges of late-type galaxies; it is therefore important to include 
late-type bulges in our analysis. Unfortunately, we do not have 
a direct measurement of the age-distribution of these bulges 
as a function of their velocity dispersion. However, bulges lie 
very close to the fundamental plane of elliptical galaxies, and 
appear to have similar formation epochs (Falcon-Barroso et al. 
2002); we therefore make the reasonable assumption that the 
age-distribution of a bulge is identical to that of an early-type 
galaxy with the same velocity dispersion. Furthermore, we fol- 
low Sheth et al. (2003) and identify a v = V c /\/2 as the bulge 
velocity dispersion for a late-type galaxy with circular veloc- 
ity V c . This identification is justified by a recent, more detailed 
population synthesis study (Tundo et al. 2006). 

We next assume that the stars in each of the nearby spheroid 
formed in a single burst, with the age t e identified as the mass- 
weighted stellar age, inferred from the spectral fits of early- 
type galaxies discussed above. Both the mean age, and the 
shape of the age-distribution, are found to be a strong function 



of the velocity dispersion dp/dt e = dp/dt e (a v ,t e ) (see Figure 1 
above). 

We also assume that the velocity dispersion a v remains con- 
stant after the formation of the galaxy. Dry mergers between 
gas-poor early type galaxies may indeed cause only a mod- 
est change to a v , as expected based on idealized models and 
confirmed by numerical simulations (e.g. Nipoti et al 2003, al- 
though the precise conclusions depend on the merging galaxies' 
structure and their orbital parameters; Robertson et al. 2006, 
Boylan-Kolchin et al. 2006). If the BH masses add linearly, 
then such dry mergers cannot maintain the Mbh - &v relation, 
implying that mergers are either infrequent, or that they are ac- 
companied by dissipation. Furthermore, in the high-resolution 
N-body simulations of Zhao et al. (2003a,b), dark matter halos 
grow in two phases. An initial rapid merger of many progeni- 
tors with comparable masses build up the potential well; this is 
followed by a slower accretion of lower-mass satellites adding 
mass primarily to the outskirts of the halo, without affecting 
the central structure. Our modeling implicitly requires only that 
most of the BH mass is assembled (possibly in multiple merg- 
ers) within a small fraction of the Hubble time. 

Under the assumption that a SMBH forms concurrently with 
the stars in each early-type galaxy, we then adopt the SMBH 
mass given by the relation, 

log (^) =a + piog (^_) (2) 
b \M Q J r V200kms-V 

with a = 8.15 and f3 = 3.83. This is the mean relation obtained 
by Tundo et al. (2006) while allowing an intrinsic scatter in 
Mbh at fixed a v (determined to be ~ 0.2 dex). Although cur- 
rent SMBH samples suffer from selection bias, Bernardi et al. 
(2006b) suggest that the M b h - cr v relation is unlikely to have 
been strongly affected. Wyithe (2006) recently found that an 
additional quadratic term in the right-hand side of equation (2), 
/32[log(er/200km s" 1 )] 2 , is required. Although the statistical 
significance of a nonzero > is modest (m 1 - 1 .5<r depend- 
ing on the sample used), we have repeated our analysis below 
using the best-fit model in Wyithe (2006), with (a, (i^fii) = 
(8.05,4.2,1.6), and have verified that our predictions for the 
quasar LF change by much less than their uncertainties. 

We assume that the SMBHs gain most of their mass in a lu- 
minous short-lived phase, at the end of their rapid assembly. 
During this phase, they radiate at some fraction /ecU of the Ed- 
dington limit of their final mass, for a duration ?q years. During 
this time, the SMBH is identified as a luminous quasar with a 
fixed luminosity Lq. We expect the lifetime to be of order the 
Salpeter time, fQ <~ 4 x 10 7 yr (Martini 2004; see also Shankar 
et al. 2004). Note that this lifetime is identified as the cumula- 
tive duration of quasar activity (i.e., if the activity is episodic, 
then it represents the sum of all episodes), and also assumes 
that the quasar is unobscured by gas or dust during the epoch 
of bright activity. This latter assumption may be justified in the 
scenario in which the lifetime is determined by feedback: the 
bright quasar epoch corresponds to an expulsion of the gas and 
dust from around the nucleus (which removes the obscuration 
on one hand, and terminates quasar activity on the other). This 
type of model of quasar activity has been proposed recently 
by Hopkins et al. (2005, and references therein), although we 
caution that in their model, the BHs spend a longer time in an 
obscured phase (see discussion in § 5 below). 

While it is useful to think in terms of separate Mbh _ c v and 
LQ-Mbh relations (the latter specified by /Edd), both of which 
have been empirically constrained in various galaxy samples, 
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our model predictions depend only on the combination of these 
two, i.e on the dependence of the quasar's luminosity Lq on a v . 
We therefore introduce the parameter aL giving the normaliza- 
tion of the Lq - cr v relation for a SMBH in a host with velocity 
dispersion 200 kms" 1 . We then have the bolometric luminosity 



:a L +12.72 + /?log 



(3) 



, 200km s 

where 10 4 57 L Q /M Q is the Eddington luminosity, so that a/, = 
corresponds to the Eddington luminosity of a 10 815 M Q black 
hole, i.e. to the adopted normalization of the M\^-a v relation 
and / E dd =1- In our analysis, we allow the normalization oll 
to vary with redshift; however, we keep the logarithmic slope 
fixed, <f logLQ/c/ logo\, = (i = 3.83. 

Under the above assumptions, and in the limit of a negligible 
scatter in Lq at fixed a v , the space density of luminous quasars 
at redshift z is predicted to be 



dN° 
c/Lq 



dN\ 
da v J 



dp , 



JLq 



(4) 



In the above equation, t e is the lookback time to redshift z, 
dp/dt e is unresolved and assumed to be a constant in each age 
bin of width dt e = lGyr, and Lq is the bolometric luminosity 
given by equation (3). 

Since the M^h - cr Y relation and f E dd are both likely to have 
intrinsic scatter, we allow a probability distribution dp/dLQ of 
quasar luminosities at fixed a v . To be specific, we assume that 
dp/dLQ follows a log-normal distribution, i.e. dp/dlogLQ is 
described by a Gaussian, 1 with a width independent of Lq. This 
can be achieved by allowing a L to have a Gaussian distribution 
with a mean (oil) and standard deviation Ac*/,. 2 Aai can there- 
fore be regarded as the combined scatter in the Ma - o\< relation 
(at fixed (T y ) and in / E dd (at fixed Mbh). To predict the quasar LF 
in the presence of scatter, the quantity in equation (4) has to be 
convolved with dp/dLQ, 

dN r , dN° dp 
dlogL Q - 



d\ogL Q J Q 



f 

Jo 



d\ogL\ 



'jlogL Q d\ogL Q 

dN" dp 
d/logLp da L 



(5) 



In these equations, logLf-, represents a mean logarithmic lu- 
minosity, computed from equation (3) with a L = (a L ), and in 
the last equation, the Gaussian dp/da L is evaluated at a L = 
logL Q -logL Q . 

In summary, the parameters of the model described above 
are tQ, (oti), and Act/,, describing the duration of quasar activ- 
ity, and the mean luminosity and its scatter during the quasar 
epoch of a SMBH in a spheroid with fixed a v . We will fit these 
three parameters simultaneously and independently in each red- 
shift bin, i.e. we do not require ab-initio that these parameters 
be redshift-independent. We note again that for comparison to 
previous work, it can be useful to regard the Lq - a v relation we 
use in our fits as the combination of an Mbh - c v relation and the 
Eddington ratio / E dd- 

4. COMPARING PREDICTED AND OBSERVED QUASAR LFS 

Our main task is to compare the predictions of the above 
model to the observed quasar LF. The bolometric quasar LF has 

1 Note that here and elsewhere in the paper, log refers to base 10 logarithm. 

2 Note that with this definition (Lq) / . Our definition is consistent with 
Tundo et al. (2006), i.e. equation (2) describes the logarithmic mean (logM bh ). 



been recently inferred from a compilation of multi-wavelength 
data (Hopkins et al. 2006). Here we adopt the observational 
data-points, with their quoted statistical errors, as reported in 
this compilation, but re-binned to correspond to the LF in 12 
different redshift intervals, centered at look-back times of t = 
1.5- 12.5 Gyr, in increments of 1 Gyr. Each redshift bin has 
a width corresponding to Af = ±0.5 Gyr. These data, together 
with their errors, are shown in both Figures 3 and 4. 

In order to compare our model predictions to these data, sev- 
eral additional uncertainties need to be taken into account. In 
determining the fraction of galaxies, at fixed a v , that fall within 
a given age bin, there are two main sources of errors. First, there 
is a Poisson sampling error 5n v , which is significant in those 
age- and cr-bins that have only a few galaxies (the number of 
galaxies is indicated by the label in each panel in Fig. 1). Sec- 
ond, errors in the age determinations (increasing from 8t e < 1 
Gyr for t e <, 7 Gyr to St e w 3 Gyr for t e ^ 10 Gyr) can shift 
galaxies between age-bins. We estimated this error, Sn t , by 
performing 20 Monte Carlo realizations of each of the age- 
distributions, generating a random age of each galaxy within its 
allowed range returned by MOPED. We maintain the prior that 
the age of any galaxy, when extrapolated to z = 0, cannot ex- 
ceed 13.7 Gyr, and identify 5n t as the range that includes 14 out 
of the 20 realizations. The resulting uncertainties 5n t are shown 
by the error-bars in Figure 1 . We add 8n p and 5n t in quadrature, 
to obtain the total uncertainty in the fraction of galaxies falling 
in each age-bin. 

In obtaining the total (normalized) space density of quasars, 
there is an additional uncertainty, from the measurement errors 
of the local velocity function of early type galaxies (Sheth et al. 
2003). The statistical errors are small ( <, 15%, although they 
increase toward the bright end), and systematic errors, due to 
including/excluding low S/N galaxies, are somewhat larger, 20- 
25%. We chose to include an additional, independent 15% er- 
ror in each luminosity bin in the quasar space density, although 
we note that in practice, these errors are much smaller than the 
combination of the other two sources of error and the observa- 
tional errors, discussed in the preceding paragraphs. 

Finally, there is additional uncertainty in the prediction of the 
quasar LF at the faint end (corresponding to the threshold a v rts 
100 km/s) in those models with a large scatter (Aol <; 0.5). 
In these models, the velocity function, and the age distribu- 
tions, have to be extrapolated to a v < 100 km/s. While the 
uncertainty introduced by this extrapolation cannot be rigor- 
ously quantified, we include it in our error budget as follows. 
We first compute the quasar LF in the absence of scatter, Nq = 
LQ(dN° /dLo), by extrapolating the power-law slope implied 
by equation (4), and convolve Nq with the given scatter Aa^ 
following equation (5). We then repeat this calculation, but re- 
ducing/increasing dp/dt e in each t e bin by its la uncertainty 
as discussed above. This results in a flatter/steeper slope for 
Nq, and yields a lower/higher Nq. We adopt these values as 
the la lower/upper limits on A^q. In practice, this increases the 
Poisson error on ANq significantly in models with scatter rel- 
ative to those with no scatter. As an example, at z = 0.31, at 
the faintest L the uncertainty is a factor of ~ 3 without scatter 
(shown in Fig. 3), but is about a factor of <~ 6 for the model with 
Acxl = 0.6 (shown in Fig. 4). 

In summary, given a set of parameters, fQ, (a L ), and Aa L , 
we compute the expected space density of quasars in each loga- 
rithmic luminosity bin, Nq as described in the previous section, 
and at each redshift and luminosity, we compute the error ANq 
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FIG. 3. — This figure shows the predicted and observed quasar LF at 12 
different redshifts, as labeled (corresponding to look-back times of 1.5—12.5 
Gyr, in increments of 1 Gyr). The solid (red) curves follow from a model 
with Lq — <t v normalization «l = (which corresponds to fen = Lq/Le(M = 1 
with our adopted My, — ov normalization). The model assumes no scatter, and 
a quasar phase that lasts for a Salpeter time at a radiative efficiency of 10 
percent (fQ = 4 X 10 7 years). The dashed (blue) curves indicate the uncertainty 
in the model prediction, arising from Poisson sampling, age-determination, 
and local galaxy abundance errors (in the top left panel, only the upper limit is 
visible). The black data-points have been adapted from a recent compilation 
by Hopkins et al. (2006). 



by adding in quadrature the Poisson sampling errors, the er- 
rors caused by the age-determinations, and the uncertainties in 
the local galaxy abundance. We then compute a x 2 f° r eacn 
data-point in the Hopkins et al. (2006) compilation that falls 
within the predicted range of luminosities. Here x 2 is defined 
as the ratio ANq^/ ANq^j, where dNQ.dift is the difference be- 
tween the predicted and measured quasar space densities, and 
flWQ.dift is the sum in quadrature of the modeling errors and 
the observational errors discussed above. Note that the num- 
ber of data-points available for the comparison of predictions 
vs. observations (and hence the number of degrees of freedom 
in the fit) varies with (ar). This is because <r v lies in the range 
100 km s" 1 <S a v ^ 280 km s , but the corresponding range of 
Lq shifts with a/,- 

At each redshift, the above procedure yields a total \ 2 , summed 
over the available observational data-points, as a function of ?q, 
{oil), and Aol- As a single figure of merit for the goodness- 
of— fit of a model, we further sum the total \ 2 over a U of the 
12 redshift bins, and obtain a cumulative reduced \ 2 (where the 
number of degrees of freedom, i.e. the number of data-points, 
minus 3 for the 3 independent fitting parameters, is typically of 
order - 100). 

5. RESULTS AND DISCUSSION 

In Figures 3 and 4, we show the predicted quasar LF at 12 
different redshifts, compared to the observed bolometric LF 
adapted from the compilation by Hopkins et al. (2006; shown 
as black points with error-bars). The twelve different panels 
correspond to 1 Gyr-wide bins in look-back times, centered on 




FIG. 4. — This figure shows the predicted and observed quasar LF at 12 
different redshifts, together with their uncertainties, as in Figure 3. The curves 
here assume oll = —0.3 (corresponding to a constant, redshift-independent 
mean Eddington ratio of (log /aid) = — 0.3 with our adopted Mbh — <r v normal- 
ization), a scatter of 0.6 dex in Lq at fixed <t„, and a quasar lifetime of 8 X 10 7 
years. This model a statistical fit with a reduced x 2 of 11 (see text for discus- 
sion), for the entire redshift range as a whole. Note that the presence of scatter 
increases the uncertainties relative to Figure 3. 



t e = 1.5 — 12.5 Gyr, in increments of 1 Gyr, and to redshifts as 
labeled. The curves in Figure 3 have been obtained in a "basic" 
model with the most naive set of assumptions, namely oll = 
with no scatter, and a quasar phase lasting for a Salpeter time 
at a radiative efficiency of 10 percent (?q = 4 x 10 7 years). This 
model is shown for illustration. The solid (red) curves show the 
predicted quasar LF, and the dashed (blue) curves bracket the 
estimated la uncertainties, as discussed in § 4 above. 

As the figure shows, this naive model under-predicts the 
overall normalization of the LF (especially at low redshifts), 
and predicts a shape that declines too steeply at the bright end 
(especially at high redshift). Therefore, under the most naive 
assumptions, the age-distributions we infer for the early-type 
populations do not "map" the local population of early-type 
galaxies onto the evolving quasar population. The age-distributions 
have the right property, in that larger early-types are older (as 
inferred by numerous previous works on "down-sizing"). This 
is in accord with the increase in the characteristic quasar lumi- 
nosity toward high z (e.g. Cattaneo & Bernardi 2003; Haiman, 
Ciotti & Ostriker 2004). However, over the range of velocity 
dispersions 100 km s" 1 < a v < 280 km s -1 for which we can de- 
termine the age-distribution, the corresponding quasar LF flat- 
tens considerably from z = 0toz = 2- 3. This flattening is not 
reproduced by the dependence of the observed age-distribution 
on cr v . 

We have searched through the full 3-dimensional parame- 
ter space (?q, {oil), Aaz.) for a set of values that corrects these 
deficiencies. Changes in ?q and amount to shifting the pre- 
dicted LF along the vertical or horizontal axes in Figure 4, re- 
spectively, without changing the predicted shape. An increased 
scatter (higher Aol) also shifts the LF up along the vertical 
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FIG. 5. — This figure shows, in each redshift bin, contours of constant re- 
duced x 2 in that bin. The blue/yellow/red contours correspond to \ 2 = 2.0, 1 .5, 
and 1.0. The x-axis in each panel shows the quasar lifetime in units of the 
Salpeter time, 4 X 10 7 yr, and the y-axis shows the mean at (equivalent to the 
mean Eddington ratio with our adopted Mbh - <x v normalization). The scatter 
in all models are fixed at Act/, = 0.9. 



axis, but it also flattens the predicted LF, by preferentially in- 
creasing the abundance at the bright, steep tail. 

The solid (red) curves in Figure 4 show predictions from a 
model with parameter choices of (fQ, (ai), Aa^) = (8 x 10 7 
years, -0.3, 0.5). This model has a reduced \ 2 of 1.1 (\ 2 = 106 
for 98 degrees of freedom over the twelve redshift bins). The 
dashed curves indicate the uncertainty in the model predictions 
in each age bin. The uncertainty increases at the faint end of 
the LF, and at lower redshifts (as well as in the highest z bin), 
as there are only a handful of SDSS galaxies in these velocity 
dispersion and age bins (see Figure 1). As noted above, the 
predicted LF at the faint end has an additional uncertainty when 
the scatter is large, since the prediction involves extrapolating 
to cr v < lOOkm/s. 

This model, while simple, is in reasonable overall agreement 
with the observational data. There does not appear to be com- 
pelling evidence that either the duty cycle or the Eddington ra- 
tio (or Mbh _ c v relation) evolves significantly with redshift. It is 
also worth pointing out that the quasar lifetime preferred by this 
model is about a factor of two longer than the fiducial Salpeter 
time of 4 x 10 7 years (at a radiative efficiency of e = 0.10). 

In Figure 5, we assume a larger scatter of Aai = 0.9, and we 
show contours of constant reduced x 2 as a function of the other 
two parameters. The reduced % 2 was computed separately in 
each of the redshift bins (each bin having ~ 10 d.o.f. on aver- 
age). The blue/yellow/red contours correspond to \ 2 = 2.0, 1.5, 
and 1 .0, respectively, and the star in each panel marks the model 
that best fits when all 12 redshift bins are combined, (?q, (a/.), 
Aa^) = (6.3 x 10 7 years, -0.5, 0.9). As this figure shows, no 
combination of parameters can match the observed LF at the 
precision of \ 2 < 1 -5 simultaneously at all redshifts (the worst 
fit is at z = 0.95, as we discussed above). On the other hand, the 
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FIG. 6. — This figure shows contours of constant reduced x 2 combined 
over all 12 redshift bins. The scatter is fixed at four different values (Aat = 
0,0.3,0.6,0.9). In each panel, the large cross marks the best-fit model. Note 
that the best-fit models have smaller reduced x 2 when the scatter is increased. 
The contour levels in each panel are adjusted to the minimum of x 2 , and are 
different in each panel. The blue/yellow/red contours enclose points with x 2 
below (4.0, 5.5, 8.0), (3.0, 5.0, 7.0), (2.0 3.0 4.0), and (1.0 1.5 2.0) in the 
four panels, respectively (from top left, clockwise). These contours serve to 
illustrate the direction of the degeneracies discussed in the text, and show that 
a non-negligible scatter (A > 0.6) is required. 



figure shows significant degeneracies: roughly, a longer life- 
time can be "traded" for a smaller oil (or equivalently, a smaller 
Eddington ratio). This is expected: a smaller oll shifts our pre- 
dicted LF curves to the left (a given Lq then corresponds to a 
larger er v , i.e. a smaller dn/da v x dp/dt e ), while a longer life- 
time moves the curves upward (scaling linearly). This explains 
the sense of the degeneracy for a given Lq. Likewise, we found 
(not shown on the figure) that an increased scatter can be com- 
pensated by a shorter lifetime and/or a smaller Eddington ratio. 

Overall, there is a suite of models that can fit the LF over 
most of the redshift range. In particular, models along lines 
with log?Q+ (ai) s» constant are degenerate, with models with 
longer lifetimes only mildly disfavored. For example, the mod- 
els with (f Q , (a L ), Aa L ) = (6.3 x 10 7 years, -0.5, 0.9) or (1.6 x 
10 8 yr, -1.6, 1.2) both have overall reduced x 2 & L0 (with all 
twelve redshift bins combined). We find that models with even 
larger scatter can have even smaller reduced % 2 ; this is a re- 
sult of the large prediction uncertainty in these models. Such 
large-scatter models, with somewhat longer lifetimes, however, 
require exceedingly small oll, or equivalently, small Eddington 
ratios ( <S 0.03), and are disfavored by other considerations (see 
discussion below). 

In Figure 6, we show contours of constant reduced \ 2 with a ll 
12 redshift bins combined. The four panels fix the scatter at four 
different values (Ac*/, = 0,0.3,0.6,0.9). In each panel, the large 
cross marks the best-fit model. As explained above, the models 
with large scatters have best-fits with lower reduced x 2 , shown 
in each panel. As the figure shows, Aa^ < 0.6 yields poor fits. 
We therefore find that the fits require a non-negligible scatter; 
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this is a result of the quasar LF being flatter at the bright end 
then the velocity-function of spheroids (as seen clearly in Fig- 
ure 2). The contour levels in each panel are adjusted to the mini- 
mum of x 2 , and are different in each panel. The blue/yellow/red 
contours enclose points with \ 2 below (4.0, 5.5, 8.0), (3.0, 5.0, 
7.0), (2.0 3.0 4.0), and (1.0 1.5 2.0) in the four panels, respec- 
tively (from top left, clockwise). These contours serve to illus- 
trate the direction of the degeneracies discussed above. 

We note that our conclusions above are based on the mass- 
weighted stellar age as a proxy for the time the SMBH forms. 
The mass-weighted age is likely to be the appropriate estimate, 
if the bulk of the stars and the SMBH forms concurrently. How- 
ever, to check the sensitivity of our results to this assumption, 
we have repeated our analysis, replacing the mass-weighted 
age either by a luminosity-weighted age, or an alternative age- 
estimate based on the (Hyp) absorption line index. In both 
cases, we find that the ages are decreased (the decrease is more 
pronounced when the absorption-line ages are used). This in- 
creases the predictions for the quasar LF at low redshifts (z 1), 
and, as a result, we find that matching the observed LF would 
require that either A«i, and/or ?q decline from z ~ 2 to 
z = by a factor of a few. One or more of the three model 
parameters would therefore have to evolve to fit the quasar LF 
over the entire range < z ^ 6. 

Our result also assumes that the luminous quasar phase is 
short, and it follows promptly the formation of the SMBH. The 
accretion history and the corresponding "lightcurve" of SMBHs 
can be more extended and complex. We first note that in the 
best-fit models we find with reduced x 2 ~ 1. the quasar lifetime 
is a factor of 1.5-2 longer than the Salpeter time, while the 
mean Eddington ratio is 0.3 - 0.5. At fixed radiative efficiency, 
this implies that during the quasar phase, the BH mass, and 
the accompanying luminosity, increases by a factor of ~ 2. In 
our models, we ignored this increase, and assumed a constant 
Lq during the luminous phase. In practice, such a variation in 
Lq would mimic a smaller Eddington ratio and a larger scatter 
in our models. Given the uncertainties about the actual quasar 
light-curve, we chose not to model this variation here. 

Observations of X-ray spectra (e.g. Wall et al. 2005; Barger 
et al. 2005) suggest that a SMBHs can suffer significant ob- 
scuration, an effect that increases at lower luminosities. For 
example, Cattaneo & Bernardi (2003) included a luminosity- 
dependent obscuration in their predictions for the quasar LF. In 
their model, the evolution of bright quasars is determined by the 
age-distribution, but the evolution of faint quasars is driven by 
the assumed form of the obscuration (which reduces the abun- 
dance of faint quasars). They showed that the obscuration can 
achieve the required flattening of the LF (which, in our case, is 
attributed to scatter instead). On the other hand, the X-ray ob- 
servations, and also theoretical arguments (Hopkins et al. 2005, 
2006) suggest that SMBHs spend a significantly longer period 
in their obscured phase, during which they shine well below 
the Eddington luminosity. Such a luminosity-dependent life- 
time would further change the shape of the predicted LF, and 
would likely require the introduction of a more detailed model, 
in which ?q, (ai) and Aaz, depend on the SMBH mass. 

Taken at face value, however, our results imply that the duty 
cycle, the Eddington ratio, and the Ma - c v relation, do not 
necessarily evolve with redshift. Our results do not, conversely, 
imply that evolution is ruled out. For example, one could adopt 
different points from the degenerate innermost red contours shown 
in Figure 3 at each redshift. Alternatively, Figure 6 shows that 



the best-fit values of (a/,) and fg do change as scatter is varied; 
one could therefore postulate an evolving (a L ) (say, decreasing 
with redshift) together with an evolving scatter (which would 
have to then increase with redshift). 

There have been suggestions in the recent literature that the 
normalization of the Ma - c v relation relation may increase to- 
ward high z, with M b h larger at fixed a v (e.g. Shields et al. 2006; 
Wandel 2004) - by as much as a factor of 4 at z = 0.36 than at 
z = (Woo et al. 2006), although these results could arise from 
a selection effect of preferentially missing the lower luminosity 
BHs at high z, or the high luminosity BHs at low z (Bernardi et 
al. 2006b). 

Eddington ratios are difficult to determine empirically (al- 
though this could be possible in the future with high precision, 
if gravity waves can be detected by LISA and identified with 
luminous quasars; Kocsis et al. 2006). Kollmeier et al. (2006) 
combined measurements of line widths and continua of a sam- 
ple of AGN at 0.3 < z < 4 to infer (/ Edd ) — 0.3, with a scatter 
of <~ 0.3dex, with no apparent evolution in the mean or the scat- 
ter (note that if the fueling rate falls below the Eddington rate, 
it would be surprising if it had no scatter). Vestergaard (2004) 
estimated Eddington ratios in a sample of high redshift quasars 
using an observed correlation between the size of the broad line 
region and the luminosity of the quasar (the correlation is cal- 
ibrated using reverberation mapping of lower redshift objects; 
e.g. Kaspi et al. 2000; Vestergaard 2002). They find values 
ranging from ss 0.1 to <; 1, with the z JSb 3.5 quasars having 
a somewhat higher mean L/L^a, and a narrower distribution, 
than the lower redshift population. In a more extensive lower 
redshift < z < 1 sample, Woo & Urry (2002) also find higher 
Eddington ratios towards z = 1 . However, these results may rep- 
resent a trend towards higher ratios at higher luminosities, and 
could also be explained by selection effects. Whether there ex- 
ists any trend, and whether it is primarily with redshift or lumi- 
nosity is an important question, but large scatter and selection 
effects presently preclude a firm answer. 

Our fits are in good agreement with these Eddington ratio 
measurements. Our results suggests that, in order to flatten the 
predicted quasar LF at z = 2 - 3, a scatter must be present in the 
relation between L and <r v . Tundo et al. (2006) find an intrinsic 
scatter of w 0.22 dex, while Wyithe (2006) and Marconi and 
Hunt (2004) both find 0.3 dex, in theM^-ay relation, and Lapi 
et al. (2006) find a similar scatter of 0.3 dex in M B h at fixed 
virial mass. The combined scatter in the M\^-a v relation (0.2- 
0.3 dex) and in /Edd (0.3 dex from Kollmeier et al.), when added 
in quadrature, is <~ 0.4 dex, which is in fair agreement with our 
finding that ^ 0.6 dex is required (for a reduced \ 2 of 1.1 or 
smaller, see Figure 6). However, it is interesting to note that our 
minimum scatter is slightly larger than the combined scatter in 
the Mbh - ov relation and /Eddd, and a factor of 3 larger than the 
scatter in the L» - sigma v relation of early-type galaxies, which 
is ~ 0.2dex. 

6. CONCLUSIONS 

In this paper, we attempted to predict the evolution of the 
optical quasar luminosity function, over the redshift range < 
z ^ 6, in a simple model with three free parameters, which is 
based on the inferred age-distribution of early type galaxies, 
and makes the assumption that quasar black hole growth and 
star-formation track one another at all times. Our "best-fit" 
models have quasar lifetimes (« 6-8 x 10 7 yr), mean Edding- 
ton ratios (w 0.3-0.5) that are in good agreement with other 
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determinations. Our results also require a non-negligible scat- 
ter between the velocity dispersion of the galaxy and the mass 
of its resident BH (rs 0.6-0.9 dex). Our basic conclusion is 
that the quasar LF can be fit to satisfactory accuracy in these 
simple models, without a compelling need for any of the model 
parameters to evolve with redshift between < z & 6. This 
result supports the direct connection between the build-up of 
spheroids and their nuclear SMBHs. 
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